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We apply the pseudoparticle approach to SU(2) Yang-Mills theory and perform a detailed study 
of the potential between two static charges for various representations. Whereas for charges in the 
fundamental representation we find a linearly rising confining potential, we clearly observe string 
breaking, when considering charges in the adjoint representation. We also demonstrate Casimir scal- 
ing and compute gluelump masses for different spin and parity. Numerical results are in qualitative 
agreement with lattice results. 

PACS numbers: 11.15.-q. 



I. INTRODUCTION 
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A common approach to gain some understanding of 
the vacuum structure of Yang-Mills theory and QCD is 
to consider effective theories restricting the path integral 
to a small subset of gauge field configurations, which are 
supposed to be of physical importance. Well known ex- 
amples are ensembles of singular gauge instantons (cf. [l| 
and references therein), which arc able to explain many 
phenomena on a qualitative level, in particular chiral 
symmetry breaking. Confinement, however, is absent in 
these ensembles, although there has been some specula- 
tion that very large instantons have the ability to cure 
this problem (0] quoted in [^). More recently related 
models for Yang-Mills theory have been proposed, which 
exhibit clear signs of confinement. There are ensembles 
of regular gauge instantons and merons 0, [3j 11] j there is 
the pseudoparticle approach @, 0, [l] , there are models 
of calorons with non-trivial holonomy [T^j and there is 
an ensemble of dyons [ll[ . The successes of these models 
regarding confinement have either been attributed to the 
long range nature of their building blocks or to maximally 
non-trivial holonomy. 



However, a satisfactory model for Yang-Mills theory 
should not only exhibit a linearly rising fundamental po- 
tential, but also also Casimir scaling for higher represen- 
tations as well as iV-ality dependence. In particular the 
adjoint potential should exhibit string breaking, when 
the corresponding charges are separated beyond a cer- 
tain distance. For a review regarding the confinement 
problem in Yang-Mills theory we refer to p^ . 



The goal of this paper is to demonstrate that the pseu- 
doparticle approach applied to SU(2) Yang-Mills the- 
ory correctly reproduces the potential between two static 
charges for various representations. Particular emphasis 
is put on string breaking in the adjoint representation. 



When computing the fundamental representation po- 
tential, Wilson loops are appropriate observables. In 
agreement with common expectation the resulting po- 
tential is linearly rising both in lattice calculations (cf. 
e.g. [l^) and in the pseudoparticle approach (cf. e.g. 0]). 
In contrast to that, the adjoint potential is expected to 
saturate at large separations, due to screening by glu- 
ons. However, various lattice studies [ij, [l^, [itI. llSl. 
[TqI . [20| . [2H have shown that the computation of this po- 
tential fails, when using Wilson loops only. One also 
obtains a linearly rising potential, because of the poor 
overlap between "string trial states" forming the Wilson 
loops and the ground state, which is essentially a "two 
gluelump state" for large separations (an exception to 
that is a rather recent study considering 3d Yang-Mills 
theory, where Wilson loops of extremely large temporal 
separation have been computed [H,!!!]). A possible way 
to overcome this problem is to consider whole sets of trial 
states including both string states (large ground state 
overlap for small separations) and two gluelump states 
(large ground state overlap for large separations). The 
adjoint potential can then be computed from the cor- 
responding correlation matrices via standard variational 
techniques. 

The paper is organized as follows. In section [TTl we 
briefly summarize the basic principle of the pseudopar- 
ticle approach, which is discussed in more detail in 0]. 
Then we compute "pure Wilson loop static potentials" 
for charges in various representations (section |III[) . We 
determine the fundamental string tension to set the scale, 
and wc demonstrate Casimir scaling for higher represen- 
tations. In section [TVl we discuss gluelump creation op- 
erators and their quantum numbers. Wc also present 
numerical results for a couple of gluelump masses and 
give estimates regarding the string breaking distance. In 
section |V] we compute the adjoint potential using both 
string-like and two-giuelump-like trial states. As ex- 
pected, it is linearly rising at intermediate separations, 
but saturates at large separations. Moreover, we perform 
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a mixing analysis showing clear evidence for string break- 
ing. Both the string breaking distance and the shape of 
the potential and its first two excitations arc in quali- 
tative agreement with lattice results. Finally we give a 
summary and a brief outlook (section IVT)) . 



II. THE PSEUDOPARTICLE APPROACH IN 
SU(2) YANG-MILLS THEORY 

A. Introduction to the pseudoparticle approach 

In the following we briefly review the pseudoparticle 
approach and its application to SU(2) Yang-Mills theory. 
For a more detailed presentation we refer to • 

The basic idea of the pseudoparticle approach is to re- 
strict the Yang-Mills path integral to those gauge field 
configurations, which can be obtained by a linear super- 
position of a small number of localized building blocks 
(pseudoparticles). A suitable choice, which is able to re- 
produce many essential features of SU(2) Yang-Mills the- 
ory, particularly a linearly rising fundamental potential, 
is given by 
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where A is the pseudoparticle size. 
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vli.v = ^6;"^ ^ hfiSoiy + Sb^So^i- The gauge field configura- 
tions entering the path integral are of the form 
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where z{i) G M.^ denotes the randomly chosen, but fixed 
position of the z-th pseudoparticle, A(i) € M its ampli- 
tude and C"^(z) S S0(3) its color orientation. Ensem- 
ble averages of observables O are defined by an integra- 
tion over pseudoparticle degrees of freedom namely am- 
plitudes and color orientations: 
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Each gauge field configuration is weighted by e ^ , where 
S is the standard Yang-Mills action 



S = ^— / d'^xF"- F°- 



pa 
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with g being the coupling constant. Such finite dimen- 
sional integrals can be computed by applying standard 
Monte Carlo techniques. 

In 0, Q it has been shown that around 400 pseu- 
doparticles H]) to ([3]) arc sufficient to reproduce many es- 
sential features of SU(2) Yang-Mills theory. In particular 
the potential between static charges in the fundamental 
representation is linear for large separations. Moreover, 
like in lattice gauge theory the scale can be set by choos- 
ing the coupling constant appropriately. Although the 
spacetime volume has been varied by a factor of ~ 16 
with the total number of pseudoparticles kept constant, 
the dimcnsionlcss ratios x^^^/'^^^^ Tcriticai/f"'^''^ are 
essentially independent of the coupling constant and their 
numerical values are in qualitative agreement with lattice 
results (a: string tension; x'- topological susceptibility; 
Tcriticai: dcconfiuement temperature). 

Although we use instanton-like building blocks ^ and 
(121), wc would like to stress that the pseudoparticle ap- 
proach is not a semiclassical model. The intention is 
rather to approximate physically relevant gauge field con- 
figurations with a small number of degrees of freedom. In 
general these configurations are not close to solutions of 
the classical Yang-Mills equations of motion. 

As has been discussed extensively in , the successful 
qualitative modeling of Yang-Mills physics seems mainly 
be related to the long range nature of the building blocks 
(P) to ([3]) and to the transversality of ^ and 



B. Numerical setup 

Unlike in previous applications of the pseudoparticle 
approach 0, [1| we consider a 4d hypercubic space- 
time region (volume: = 5.0'*) with periodic boundary 
conditions. This allows to fully exploit translational in- 
variance and to adopt certain smearing techniques from 
lattice gauge theory (cf. appendix I A 1 1 and I A 2p . Periodic 
versions of the pseudoparticles ([T]) to Q arc obtained by 
applying a blending technique introduced in 0] applied 
to all four spacetime directions (width of the blending 
region: 0.25 x L = 1.25). 

We use 625 pseudoparticles (P) to ([3]) with A — 0.5. 
This amounts to a pseudoparticle size, which is of the 
same order of magnitude as their average nearest neigh- 
bor distance. The ratio of pseudoparticles is chosen ac- 
cording to Mnstanton ■ -^antiinstanton ■ -^akyron =3:3:2. 

We consider values of the coupling constant g between 
9.5 and 18.5. This corresponds to spacetime extensions in 
the range of 1.55 fm < L < 2.31 fm in physical units, where 
the scale has been set by identifying the fundamental 
string tension with 4.2/fm^. The majority of computa- 
tions has been performed at g ~ 12.5, which amounts to 
L = 1.85 fm. 

All ensemble averages have been computed from 50 
thermalized gauge field configurations ^ with indepen- 
dently chosen pseudoparticle positions. For efficiency of 
the computation we convert these continuum gauge field 
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configurations to lattice link configurations before ensem- 
ble averages are computed (cf. appendix 1^ . 



III. PURE WILSON LOOP STATIC 
POTENTIALS FOR DIFFERENT 
REPRESENTATIONS 

In this section we compute the potential between two 
static charges for different representations from Wilson 
loops only. This amounts to considering temporal corre- 
lations between string trial states 



(^(•^)(x))t;7('')(x;y)0('^)(y)|l]), 



(8) 



where (0('^')t and 0'^-^' represent static charges in rep- 
resentation J = 1/2, 1,3/2, ... at X and y and f/^"^' is 
a straight parallel transporter connecting the charges in 
a gauge invariant way. We orient fj^"^-* along one of the 
four space diagonals allowing us to consider rather large 
spatial separations R = |x — y| without being affected by 
periodicity. To maximize the ground state overlap, t/^'^-' 
is approximated by a product of APE smeared spatial 
links (cf. appendix lA To compute effective masses. 



as a function i 




we use the variational method explained in appendix IB] 
We consider correlation matrices built from three string 
trial states ([8]), which differ in their APE smearing pa- 
rameters (A^APE S {5 , 15 , 35}, QfAPE = 0.5). 

To give the reader an idea of the plateaux quality ob- 
tained from our pseudoparticlc computations, we show 
ground state effective masses for the fundamental rep- 
resentation potential in FIG. [1^. As potential values we 
take weighted averages of effective masses in ranges where 
plateaus are indicated (the solid horizontal lines in 
FIG. [TJi). The corresponding potential is plotted in 
FIG. [TJd. As expected it is linear for large separations. 
For separations smaller than the pseudoparticlc size and 
average nearest neighbor distance cutoff effects are dom- 
inant yielding a parabolic rather than a Coulomb-like 
behavior. To set the scale, we perform a minimizing 
fit with V{R) = Vq + aR to the data points indicated in 
FIG. [T]d and identify the resulting string tension 
cr = 0.57(3) with the "physical value" 4.2/fm^ We ob- 
tain a spacctime extension L = 1.85 fm. 

It is well known that Wilson loops in higher representa- 
tions > 1> 3/2, 2, . . ., can be expressed in terms 
of fundamental representation Wilson loops 
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FIG. 1: a) Effective masses for the static potential in tiie fun- 
damental representation for different separations R as func- 
tions of T. b) The corresponding static potential l/(^/-^' as a 
function of the separation R. 



Results for the corresponding potentials are shown in 

FIG.IIl. In FIG.Hd we plot ratios ^("'V^^^^^^- For inter- 
mediate distances these ratios arc expected to be close to 
the ratios of the corresponding Casimir operators, which 
are given by 8/3 (J = 1), 5 (J = 3/2), 8 (J = 2) and 
35/3 (J = 5/2) (the sohd hues in FIG.[2}3). While the ad- 
joint potential is in excellent agreement with the Casimir 
scaling hypothesis, higher representations exhibit certain 
deviations at larger separations. These findings are in 
agreement with both Q and [2^, where similar analyses 
in related meron and regular gauge instanton ensembles 
as well as in 4d SU(2) lattice gauge theory have been 
performed. 

Note that there is no sign of string breaking in the 
adjoint potential (J = 1 curve in FIG. [5^), which is 
expected to happen at distances R w 1.0 fm . . . 1.25 fm 
[20, [2al ■ As we will demonstrate in section |V] this is due 
to the poor ground state overlap of string trial states ([S]) 
for separations larger than the string breaking distance. 
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a) V*''' as a function of R for different J 
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b) V*-^' / V'^'^' as a function of R for different J 
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FIG. 2: a) Pure Wilson loop static potentials 1/'''' for difler- 
ent representations J as functions of the separation R. b) Ra- 
tios of different representation static potentials /'(/(i/^) 
compared to the Casimir scaling expectation as functions of 
the separation R. 



IV. GLUELUMP MASSES FOR DIFFERENT 
SPIN AND PARITY 

To estimate the string breaking distance, we compute 
masses of states containing a single static adjoint charge 
surrounded by gluons, so called gluelumps. 



A. Gluelump trial states 

The symmetry group of states constrained by a single 
static charge (j)^^^ located at position x is SO{3) x P 
(rotation and parity both with respect to x). Suitable 
"gluelump creation operators" with well defined quantum 
numbers are given by 



^(,/_i,P_+)^^^ = Tr($(x)B,(x)) 
g1'^=^-^=-)(x) 



G 



(j=2,p=-; 



(x) 



Tr(^$(x)e,Mi?fcS,(x) 
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where j e {x,y, z}, ^ = ((?!)(^^)°cr°/2, Bj denotes the 
color magnetic field and Dj = dj — . . .] the covari- 
ant derivative. In the literature ( J = 1, P = +) and 
(J = 1, P = — ) gluelumps are also referred to as mag- 
netic and electric gluelumps respectively. 

Since we perform a latticization of our pseudoparti- 
cle gauge field configurations, it is convenient to replace 
Bj by the magnetic clover leaf in a plane perpendicular 
to the j-direction and DkBi by the corresponding electric 
clover leaf in a plane perpendicular to the Z-direction [l^ . 
After this replacement our gluelump creation operators 
are not rotationally symmetric anymore, but belong to 
one of the irreducible representations of the cubic rota- 
tion group Oh: Ti for and (|14p . i.e. the correspond- 
ing states are superpositions of spin J = 1,3,4,..., and 
T2 for p5|) . i.e. spin values J = 2, 3, 4, . . . are possible. 



B. Numerical results 

We determine gluelump masses from temporal corre- 
lations 



(O|(Gf^)(x,r))^Gf^)(x,0)|l^) ( 



no sum over j ) , 

(16) 

where the spatial links appearing in the clover leafs are 
APE smeared, to maximize the ground state overlap 
(TVape = 5, aAPE = 0.5; cf. appendix I A 1 1 and [Bjl . 

Note that gluelump masses by themselves are not phys- 
ically meaningful. They are cutoff dependent quantities, 
which diverge in the continuum lirnit, due to the self 
energy of the static charge (cf. e.g. [H, [131 )■ However, 
mass differences of gluelumps or more generally of states 
containing the same number of static charges arc phys- 
ical obscrvables. This offers the possibility to compute 

mass differences between gluelumps and to estimate the 

(J P) 

string breaking distance Rl^^" ' by intersecting the pure 
Wilson loop adjoint potential F^^-* with two times the 
corresponding gluelump mass (cf. FIG. [3^). Results are 
collected in TABLE HI 

Comparing these results to available lattice results we 
find rather good agreement for the estimated "magnetic 
string breaking distance": ^gbil^ltico = 1.0 fm ... 1.25 fm 
[20I. [25!. Moreover, the mass of the magnetic and of the 
electric gluelump are of the same order of magnitude. 



UP) 


gluelump 


(j,p) (1,-) 

gluelump gluelump 




(1,+) 
(1,-) 

(2,-) 


1037(21) MeV 
865(17) MeV 
1123(23) MeV 


172(38) MeV 
257(40) MeV 


1.01(3) fm 
0.85(2) fm 
1.09(3) fm 



TABLE I; Gluelump masses, mass differences and estimated 
string breaking distances for different quantum numbers 
(J,P). 
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a) V* ' as a function of R 
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FIG. 3: a) Estimating string breaking distances R^^^'^^ by 
intersecting the static adjoint potential y'^' with 2 m^J^^^^_^_^^. 

b) Estimated string breaking distances fl^'j^'^^ as functions 
of the coupling constant g for different gluelump quantum 
numbers ( J, P). 



However, their mass difference, which is hard to measure 
precisely, due to large absolute mass values, has the oppo- 
site sign: (m^};+|^^p - m^l^^Lip) lattice = -203(76) McV 
[2^. An explanation might be that ultraviolet fluctua- 
tions, which are probably more important for localized 
objects like gluelumps than for the static potential at 
large separations, arc not adequately described by our 
pseudoparticle regularization (the resolution of our gauge 
field configurations is roughly given by the pseudoparti- 
cle size and the average nearest neighbor distance, which 
is « 0.4 fm for g = 12.5). Similar problems have been en- 
countered in a study of the gluelump spectrum in related 
meron and regular gauge instanton ensembles Q. 

We have also investigated the stability of the estimated 
string breaking distances under a variation of the cou- 
pling constant g (cf. FIG.[3}3). While the spacetime vol- 
ume in physical units has been increased by a factor of 
w 5 {g = 9.5 . . . 18.5 corresponds to 

L ~ 1.55 fm . . . 2.31 fm) with the total number of pscu- 
doparticles kept constant, the estimated string breaking 
distances Ril,'~^\ Rih~^ and r'"^'~^ vary by only w 6%, 



« 9% and 19% respectively. This complements the 
scaling analysis performed in [8|, where the dimension- 
less ratios x^/^/ct^/^ and Tciiticai/o'^''^ have been found 
to be essentially independent of the coupling constant. 



V. ADJOINT STRING BREAKING 

In this section we supplement our basis of string trial 
states (|8]) used in the variational method (cf. appendixlBj) 
by a second type of trial state, which resembles a two- 
gluelump state. 



A. Two-gluelump trial states 

When computing the adjoint potential, it is essential to 
have a basis of trial states with significant ground state 
overlap for arbitrary separations of the static charges. 
For intermediate separations the ground state is expected 
to be a gluonic string connecting the charges, i.e. a state 
with large overlap to the string trial states ([8]). For large 
separations the charges are supposed to be screened by 
gluons forming two essentially non-interacting gluelumps. 
As we will demonstrate at the end of this section, the 
overlap to string trial states is rather poor. 

In the following we consider the static charges ((/''■^•')^ 
and (j)^^'' located at positions x = (0,0, —z/2) and 
y = (0,0,-|-z/2). The corresponding symmetry group is 
S0(2) X P^{xP^). The S0(2) rotation is around the z- 
axis with angular momentum as corresponding quan- 
tum number. denotes reflection along the z-axis with 
respect to its center, and in the case oi Jz = there is 
another symmetry P^, reflection along the a;- axis [28j . 

The quantum numbers of the string trial states ([8]) 
are Jz = 0, Pz — + and = +. Since we are not 
only interested, whether the adjoint potential saturates 
at two times the gluelump mass, but also whether the 
string actually breaks, when static charges are separated 
adiabatically, we need two-gluelump trial states, which 
have the same quantum numbers. Using products of two 
single-gluelump creation operators (fT5|) to (fT5|) an obvi- 
ous choice is 

(G(^^^)(x)GW^)(y) + G('^-^)(x)G«^)(y) + 

GW^)(x)Gi^-^'(y))|f^)- (17) 

(J p) (J p) 

There is another possibility, where Gz Gz is 
weighted by a factor of —2. However, since string break- 
ing is closely related to non-vanishing off-diagonal ele- 
ments in correlation matrices, we expect such a state to 
be less suited for our purposes, due to the relative minus 
sign. Therefore, we supplement our basis of string trial 
states by (|17l) for the following computations. 
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B. Numerical results 

Although the electric gluelump turned out to be lighter 
than the magnetic gluelump, we compute string breaking 
both with electric and with magnetic trial states. This 
allows a direct comparison to lattice results, since there 
seem to be only investigations of 4d adjoint string break- 
ing with magnetic trial states in the literature [20, Hlj . 

We extract the adjoint potential from 3x3 correlation 
matrices containing two string trial states 
{Nape & {15 , 35}, oape = 0.5; cf. appendix lA 1|) and 
one two-gluelump-trial state (A^ape = 5, oape = 0.5) 
by means of the variational technique explained in ap- 
pendix |B] Results, which are shown in FIG. HI arc qual- 
itatively the same both for magnetic and electric trial 
states. In contrast to the pure Wilson loop computation 
(cf. FIG. [2^, J = 1 curve) the potential saturates at two 
times the gluelump mass and at separations close to the 
estimated string breaking distance (cf. TABLE [l|. Wc 
also plot the first and second excitation. It is interesting 
to note that for small separations the first excitation is 
a string excitation, for intermediate distances it becomes 
a two gluelump state and for large separations it is a 



string state again. This level ordering is in ag reement 
with various lattice computations in 3d [l^, and 4d 
[iOl SU(2) Yang-Mills theory. Moreover, within « 20% 
there is agreement with [20j regarding the string breaking 
distance and the separation of the energy levels. 

Wc have also investigated, whether the string actu- 
ally breaks, when static charges arc separated adiabat- 
ically, or whether there is just a plain level crossing of 
string and two-gluelump ground states. Since this ques- 
tion is rather hard to resolve from potential plots like 
FIG. m we perform a mixing analysis in close analogy to 
[20l |. To keep things as simple as possible, we consider 
two normalized trial states, a string trial state jstring) 
(Aape — 35, aAPE — 0.5) and a two gluelump-trial state 
two-gluelump) (A'ape = 5, oape = 0.5). From the vari- 
ational method we obtain approximations of the ground 
and the first excited state: 

|0) « a°t,i„g I string) +a°„„_gi^„i^„p I two-gluelump) 

(18) 

1 1 ) ~ flstring I String) + a J„o-giuciump I two-gluclump) . 

(19) 



a) magnetic two-gluelump trial states 



> 

0) 



3000 
2500 
2000 
1500 
1000 
500 




second excited state 
first excited state 
ground state 



3000 
2500 

> 

J 2000 
~~ 1500 
^ 1000 
500 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 
R in fm 

b) electric two-gluelump trial states 
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a) overlaps of the ground state approximation 
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b) overlaps of the first excited state approximation 
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FIG. 4: The static adjoint potential 1/'^^ and its first two 
excitations as functions of the separation R. a) Magnetic 
two-giuelump trial state, b) Electric two-gluelump trial state. 



FIG. 5: a) Overlaps of the ground state approximation to 
the trial states as functions of the separation R. b) Overlaps 
of the first excited state approximation to the trial states as 
functions of the separation 7?. 
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In FIG. [5] we show the squared amphtudes |agjj.;jjgp and 

I'^two-giuciumpl^ ^ functions of the separations for the case 
of a magnetic two-gluehimp trial state. It is clearly visi- 
ble that for separations smaller than the estimated string 
breaking distance Rsli'~^^ ~ 1.0 fm the ground state is es- 
sentially a string state, while the first excited state is 
almost exclusively a two-gluelump state. For separations 
larger than the string breaking distance the situation is 
reversed, explaining why a computation of string break- 
ing from Wilson loops alone has failed. In a narrow range 
around the string breaking distance we observe mixing of 
string and two gluelump trial states. This mixing indi- 
cates that there is a smooth transition from a string state 
for R<R^l^''^^ to a two gluelump state for R>r'^1^^\ 
when static charges are separated adiabatically. We con- 
clude that string breaking is present in the pseudoparticle 
approach. Comparing the overlap plots from FIG. [5] to 
the lattice result in [20] we find again rather good agree- 
ment. 



VI. SUMMARY 

We have performed a detailed study of adjoint string 
breaking in the pseudoparticle approach. In agreement 
with lattice gauge theory the static potential saturates 
at two times the gluelump mass, which corresponds for 
magnetic trial states to a charge separation of « 1.0 fm. 
Moreover, from a mixing analysis we have obtained 
strong indications that the connecting gluonic string ac- 
tually breaks, when the corresponding charges are sepa- 
rated adiabatically. We have also computed static poten- 
tials for various representations from Wilson loops only. 
There is excellent agreement with the Casimir scaling 
hypothesis for the adjoint potential and also higher rep- 
resentations exhibit only minor deviations. In view of 
these successes we conclude that the pseudoparticle ap- 
proach is a model for SU(2) Yang-Mills theory, which 
correctly reproduces many essential features connected 
to confinement. 

Gluelump masses, on the other hand, differ to some 
extent from lattice results. In particular the level or- 
dering of the magnetic and the electric gluelump is re- 
versed. This might indicate the limitations of pseudopar- 
ticle regularizations, where degrees of freedom are chosen 
to model long range correlations rather than ultraviolet 
fluctuations. 



Outlook 

After this successful computation of screening of ad- 
joint charges within the pseudoparticle approach, a nat- 
ural next step is an investigation of string breaking in 
QCD (for a recent lattice study cf. HI]). First steps 
regarding the treatment of fermionic fields in terms of 
pseudoparticles have already been taken [sO, [3l| . An in- 



teresting feature of the pseudoparticle approach regard- 
ing numerical efficiency is the fact that it uses a rather 
small number of degrees of freedom compared to lat- 
tice gauge theory. This makes exact computations of 
fermionic all-to-all propagators not only possible, but 
also extremely cheap. This might offer the possibility to 
compute fermionic quantities on a qualitative level with- 
out using high performance computer resources. 
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APPENDIX A: LATTICIZATION OF 
PSEUDOPARTICLE GAUGE FIELD 
CONFIGURATIONS 

Before computing ensemble averages we convert the 
continuum gauge field configurations (|4|) to 24^ lattice 
link configurations. We do this to increase the efficiency 
of our computations and to adopt certain smearing tech- 
niques from lattice gauge theory. The corresponding 
links are computed by sampling (|4]) sufRciently often 
along these links and by multiplying corresponding SU(2) 
matrices. We would like to stress that, although ensemble 
averages are computed on such lattices, the underlying 
gauge field configurations are still continuum gauge field 
configurations. Therefore, all results presented in this pa- 
per are results from a continuum model, where the lattice 
has only been introduced for the sake of convenience and 
numerical efficiency. 



1. APE smearing of spatial links 

The ground state overlaps of string, gluelump and two- 
gluelump trial states (eqns. ^ to ^ and dTT])) 
can be increased by giving the corresponding creation 
operators certain volume extensions. Such operators can 
be obtained by replacing all spatial links Uj = C/j°^ by 
their APE smeared versions [/j^*^'^^: 

Ur'\-) - PsU(2)(c/r(.)+«APE E 

\ fc=±l,±2,±3 

Ut''\x)uf\x + ek)u[1\x + ek + e,)) , (Al) 

where Psu(2) denotes an appropriate normalization pro- 
jecting back to SU(2) (3^ . 
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2. HYP smearing of temporal links 

To reduce the self energy of static charges, which in 
turn significantly improves the signal-to-noise ratio, we 
use HYP smearing of temporal links [1^ . HYP smearing 
is reminiscent to three iterations of APE smearing, where 
links outside a hypercube around the original link are ig- 
nored. There arc three parameters, which have been "op- 
timized" in [131 and are commonly referred to as HYP 2: 
t?HYP2 = (1.0, 1.0, 0.5). 

Throughout this paper we use three iterations of HYP2 
smearing for temporal links. Since the extension of the 
resulting links is still below the pscudoparticlc cutoff, 
which is of the order of the pscudoparticlc size and the 
nearest neighbor distance, we do not expect to alter cor- 
relation functions at separations, where physically mean- 
ingful results can be extracted. We checked that this 
is indeed the case for various observables. An example, 
the static potential in the fundamental representation, 
is shown in FIG. [6l It is obvious that the slope of the 
potential for large separations is the same for unsmcared 
and for HYP2 smeared temporal links, while statistical 
errors are significantly reduced for the latter. 

V*'"^* as a function of R 
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FIG. 6: The static potential in the fundamental represen- 
tation as a function of the separation R (unsmeared 
temporal links, one level of HYP2 and three levels of HYP2). 



APPENDIX B: CORRELATION MATRICES AND 
EFFECTIVE MASSES 



To determine gluclump masses and static potentials 
more reliably, we use a well known variational technique 
(cf. e.g. IH). 

The starting point is a correlation matrix 



Cjk{T) 



moj{T)] OKiom), 



(Bl) 



where Oj are suitable creation operators yielding a trial 
state basis of string, gluclump or two gluelump states 
with appropriate quantum numbers and possibly differ- 
ent extensions. The basis should be chosen such that a 
good approximation of the ground state is possible. 

Once this correlation matrix has been computed, one 
has to solve the generalized eigenvalue problem 



CMTo)v'j^HTo) = CjK{To~a)v'j^\To)X^''^ (B2) 



at a fixed value of Tq. Approximations of low lying states 
\N) are given by 



(B3) 



and corresponding energies can be determined from ef- 
fective mass plateaus: 



(N) 

effective 



(T) 



= In 



iv^f'\n)rcjKiTyK'>{To 



{vf\T,)YCjK{T-ay^'{T,) 



{N)i 



(B4) 
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